Entanglement distillation for continuous- variables under a thermal environment: 

Effectiveness of a non-Gaussian operation 



m 

o 

(N 



Jaehak Lee and Hyunchul Nha 

Department of Physics, Texas A & M University at Qatar, P. O. Box 23874, Doha, Qatar 



We study the task of distilling entanglement by a coherent superposition operation ta+rci' applied 
to a continuous-variable state under a thermal noise. In particular, we compare the performances 
of two different strategies, i.e., the non-Gaussian operation ta + ra^ is applied before or after the 
noisy Gaussian channel. This is closely related to a fundamental problem of whether Gaussian or 
non-Gaussian entanglement can be more robust under a noisy channel and also provides a useful 
insight into the practical implementation of entanglement distribution for a long-distance quantum 
communication. We specifically look into two entanglement characteristics, the logarithmic nega- 
tivity as a measure of entanglement and the teleportation fidelity as a usefulness of entanglement, 
for each distilled state. We find that the non-Gaussian operation after (before) the thermal noise 
becomes more effective in the low (high) temperature regime. 

PACS numbers: 03.67.Mn, 03.65.Yz, 42.50.Dv 
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I. INTRODUCTION 

In quantum information processing, it is an important 
task to distribute entanglement between distant parties, 
but entanglement can be easily degraded due to inter- 
action vifith a noisy environment. Numerous schemes to 
overcome the decoherence by using nondeterministic lo- 
cal operations were proposed to distill entanglement for 
discrete variable systems [ll-Q- In the continuous vari- 
able (CV) regime, it is known that Gaussian states can- 
not be distilled by using only Gaussian operations [1-Q. 
A frequently used CV entangled resource, i.e. two-mode 
squeezed vacuum (TMSV), is a Gaussian state and a cer- 
tain non-Gaussian operation is thus required. Distillation 
of Gaussian states was initially studied with two elemen- 
tary non-Gaussian operations, namely, a single photon 
subtraction a and addition a^ [7- 10]. The entanglement 
distillation by the photon subtraction scheme was ex- 
perimentally realized [ll|. Recently, it was also found 
that a more efficient distillation can be achieved by a co- 
herent superposition of photon subtraction and addition, 
ta + ra^ |12l |. 

These studies for CVs, however, did not take into ac- 
count the interaction with a noisy environment that ac- 
tually motivated the topic of entanglement distillation. 
They have simply shown that the non-Gaussian oper- 
ations can distill a pure Gaussian entangled state into 
a pure non-Gaussian entangled state with higher entan- 
glement. The task of entanglement distillation is more 
important when the state becomes mixed due to interac- 
tion with a noisy environment. Only a few works so far 
treated this practically important problem, which found 
that distillation by a single photon subtraction is possible 
even in the presence of noise [ij] . 

In this paper, we investigate how efficiently we can dis- 
till CV entanglement under a thermal environment by a 
coherently superposed operation ta + ra' . While the pre- 
vious analysis for the distillation of mixed CV entangled 
state considered the case of applying the operation after 



entangled states are distributed through noisy channels 
[ia |. we also study another strategy. That is, a non- 
Gaussian operations is applied before a CV entangled 
state is distributed to distant parties. Which of the two 
strategies gives a better performance is closely related 
to a fundamental problem of to what extent Gaussian 
and non-Gaussian entangled states can be robust un- 
der a noisy channel. Some evidences were put forward 
to support the conjecture that Gaussian entanglement 
is more robust than non-Gaussian entanglement |14l.ll5l[. 
However, there exist some counterexamples in which non- 
Gaussian states can be more robust than Gaussian states 
[16,-18] . From a practical point of view, it is crucial to 
have a longer survival time of entanglement, as entangle- 
ment cannot be distilled at all once it dies out. 

We consider a TMSV (a prototype of CV entangled 
state) as an initial state and apply ta + rd^ (an elemen- 
tary non-Gaussian operation that includes the photon 
subtraction and the addition as special cases) before or 
after the noisy channel [Fig. [Tj. We will show that two 
different strategies have advantages in different temper- 
ature regimes. A non-Gaussian entangled state distilled 
before the thermal noise can survive longer than a Gaus- 
sian entangled state without any operation. This ef- 
fect is particularly remarkable in the high-temperature 
regime where the survival time of Gaussian entangle- 
ment is short. A Gaussian state that becomes separable 
via a noisy channel cannot be distilled into an entangled 
state because local operations and classical communica- 
tion (LOCC) cannot create any entanglement from a sep- 
arable state. However, if the entanglement survives in a 
Gaussian state, it can be distilled into a highly entangled 
state by the coherent superposition operation. In the low- 
temperature regime where Gaussian entanglement sur- 
vives long enough, it turns out that the non-Gaussian 
operation after the noisy channel enhances entanglement 
better. 

This paper is organized as follows. In Section [Hi we 
briefly introduce the description of two-mode entangled 



states and their tranformations by a coherent superposi- 
tion operation and under a thermal noisy channeL Then, 
we study the entanglement properties of output states, 
namely, the logarithmic negativity in Section [IIII and the 
teleportation fidelity in Section ITVl We summarize our re- 
sults and discuss the applicability to practical protocols 
in Section FVl 



II. CV ENTANGLED STATES AND NOISY 
CHANNEL 

A. Two-mode squeezed vacuum and coherent 
superposition operation 

We start with a Gaussian entangled state generated by 
two-mode squeezing, which is described by 



scribed by a master equation 



|^)tmsv = sechs y tanh" s\n, n). 



n=0 



(1) 



where s is squeezing parameter and \n,m) represents a 
two-mode state in Fock state basis. The corresponding 
characteristic function is given by 
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where Di{a) = e"°»~" "• denotes the displacement op- 
erator acting on mode i = 1,2 with amplitude a [1^ . 

A coherently superposed operation ta + ra^ can be im- 
plemented by erasing the which-path information of the 
trigger photon emerging from the photon subtraction a 
or the photon addition a^ [20|. The parameters t and r 
can be controlled by adjusting the transmitivity of the 
beamsplitter that erases the which-path information be- 
fore photodctection. When this operation is applied to 
the mode i of a multi-mode state p, whose original char- 
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Here T and rith are the loss coefficient and the aver- 
age photon number in thermal environment, respectively. 
With current technology, the loss coefficient in optical 
fiber can be made less than a few dB/km. L[A] is the 
Lindblad operator defined by L[A]p= 2ApA^ — A^Ap — 
pA^ A. Here, we assume that both of the two modes prop- 
agate under the thermal channels with the same F and 

nth- 

Instead of solving the master equation, the interaction 
with the thermal noise can be described by a beamsplitter 
model. The input state p and the thermal ancilla state 

Pth = J _^atai "^ith mean photon number nth = j^^ 

are mixed at a beam splitter with transmissivity rj = e~ 
and the ancilla mode is then traced out in the output. 
That is, 



p' = Trth \uBsP®PthUl 



BS 



(5) 



where Trth is the partial trace over the ancilla mode and 
^BS = e^'^"i°"'~°'°tii) describes the beamsplitter interac- 
tion of ith mode and ancilla mode (cos = y/rj) . Then, 
the characteristic function for the output state turns out 
to be 

x'{'t) = N.ixCt)) 



= Xth( Vl - V^t) X 

x{£.i,--- ,^i^i,Vv£.i,^t+i,---), (6) 

where Xth(C) = e~2(2"th+i)? jg ^^q characteristic func- 
tion of a thermal state with mean photon number nth- 

Combining Eqs. (3) and (6), one can obtain the char- 
acteristic functions of final states via two different ways 
of distillation [Fig. [T] using the characteristic function of 
initial state xtmsv with operators Oi and Ni [i ~ 1,2). 
If we perform distillation before sending the initial state 
through noisy channels [Fig. [1] (a)] , the output state will 
be given by 

Xno(Cl,6) = N^ O A^2(Ol O 02(XTMSv(ei,6)))- (7) 

On the other hand, if we perform distillation after send- 
ing the state through noisy channels [Fig. [1] (b)], the 
output state will be given by 

Xo„(a,6) = Oi°02(Afio7V2(xTMSv(Ci,6)))- (8) 

In the following, we compare the entanglement proper- 
ties under two different strategies with the corresponding 
characteristic functions Xno and Xon- 



B. Noisy Gaussian channel 



III. LOGARITHMIC NEGATIVITY 



On the other hand, when the ith mode propagates 
through a thermal-noise channel, its evolution can be de- 



In this section, we first quantify the degree of entan- 
glement for each output state under the two different 



operation, we have 




FIG. 1. Schematic diagram for CV entanglement distillation 
under noisy channels 



schemes [Fig. [T]. As a measure of entanglement, we 
particularly adopt the logarithmic negativity, an easily 
computable entanglement monotone [2l|, ^^ . 



A. Estimation of Logarithmic Negativity 



The logarithmic negativity is defined as 



ENip) 



log2||p^-*||l 



(9) 



where p'^-* is the partial transpose of the density matrix 
p and the trace norm is defined as \\A\\i = try/A^A. 

For a Gaussian state, which can be completely de- 
scribed with its first and second momenta of canonical 



operators Xi 



--{a^ + al), the 



-l^{a, + a]) and p, = -^y 
logarithmic negativity can be easily calculated from its 
second-moment covariance matrix a. The logarithmic 
negativity of a Gaussian state is given by [2l| 



EN{cr) = max |o, - log2(2rf_)| 



(10) 



where d- represents the least symplectic eigenvalue of 
partially transposed covariance matrix a'^^ . 

On the other hand, the distillation operation we con- 
sider is non-Gaussian, thus Eq. (|10p is not adequate to 
fully address the logarithmic negativity by incorporating 
all higher-order momenta. For a non-Gaussian state, we 
investigate the density matrix elements in the number- 
state basis. For a numerical calculation of negativity, we 
restrict our consideration to the subspace of finite pho- 
ton number states. That is, we calculate the negativity 
with a truncated density matrix ptrunc = PopPo where 
pQ = X^n '771=0 !"■' "^)("'' "^1 '^ith the truncation number 
A'trunc- Note that the logarithmic negativity of ptrunc 
cannot be larger than that of p due to the monotonicity 
of E]s[. The trace norm Hp"^-*!]! is a monotone under a 
positive partial transpose preserving operation (PPT op- 
eration) which maps p into pi = ^i(p)/tr^i(p) with its 
probability pi = tr*j(p) ^. That is, 
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(12) 



(13) 



Since the set of operations {Pq, Pi — I — Pq} is a PPT 



This implies that the calculation with the truncated, un- 
normalized, density matrix does not overestimate the log- 
arithmic negativity of the original state. In particular, 
if p is separable, the negativity of ptrunc necessarily be- 
comes 0. Furthermore, the difference between the actual 
negativity and the truncated negativity would be negli- 
gible if the truncation number A'trunc is taken sufficiently 
large compared to the mean photon number of a given 
state. 

The density matrix elements in the number-state basis 
can be calculated using the characteristic functions [I^l 
of Eqs. dZl) and dH) as 

PtjM = {i,j\p\k,l) 

(14) 

In the following, we present numerical results using three 
squeezing parameters s = 0.029, 0.114 and 0.403 for the 
initial TMSV under different schemes. We note that in 
the pulsed-regime of squeezed light, the squeezing level 
s — 0.403 (3.5dB) is currently achievable using an optical 
parametric amplifier 23] . We consider a thermal photon 
number ?ith = 10^^ for a low-temperature reservoir and 
nth — 10^^ for a high-temperature reservoir as examples. 
When we apply a coherent operation ta + rd^ to each 
state, the value of t is numerically optimized to yield the 
highest degree of entanglement under a given situation. 
We have checked that the local mean photon number is 
less than 1 for all considered states thus restricted the 
truncation number to A^trunc = 5. Furthermore, an in- 
crease of A'trunc to a higher number did not show any 
appreciable changes in our results. 



B. Result 

In Fig. [21 we plot the logarithmic negativity as a func- 
tion of 77 = e~'"* that represents the interaction time 
with a thermal reservoir [Eqs. (4) and (5)]. When rj = 1, 
i.e. no reservoir-interaction, it is known that the coher- 
ent operation (thick curves) can enhance the entangle- 
ment more effectively than the mere photon subtraction 
(thin curves), which is remarkable particularly in the low- 
squeezing regime [l2|. This is also true for any value of 
ry < 1 as shown in Fig. [21 As the environmental inter- 
action becomes longer, i.e. ry decreases below 1, we see 
that the two strategies applying the coherent operation 
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FIG. 2. Logarithmic negativity as a function of r; = e"""', 
which characterizes the interaction time with a thermal envi- 
ronment [Eqs. (4) and (5)] under different strategies: no dis- 
tillation (black dotted), subtraction before noisy channel (thin 
blue solid), subtraction after noisy channel (thin red dashed), 
coherent operation before noisy channel (thick blue solid) and 
coherent operation after noisy channel (thick red dashed). 
Left (Right) panels represent a low (high)-temperature reser- 
voir with nth = 10~^ (10~^). Squeezing parameters are given 
by (a), (b) s = 0.029, (c), (d) s = 0.114 and (e), (f) s = 0.403. 
Two thin curves that represent the photon subtraction before 
and after the noisy channel are overlapped in (a), (c), and (e). 



before and after the noisy channel, respectively, provide 
an advantage in different temperature regimes. 

When the thermal photon number of the reservoir is 
very small (left panels), the coherent operation after the 
noisy channel (dashed curves) generally makes the output 
entanglement higher. Even for a very low rj, the output 
entanglement can maintain a rather high value by the 
coherent operation, particularly for an initially weakly 
squeezed state [Fig. [2] (a)]. Note that the value of the 
parameter t in the coherent operation ta + ra^ is opti- 
mized for each state case by case in order to maximize 
the output entanglement in all plots. 

The coherent operation or the photon subtraction 
transforms a Gaussian state to a non-Gaussian state. 
Therefore, our result implies that keeping Gaussianity 
through the noisy channel and applying later the coher- 
ent operation is advantageous for an optimal distillation 
in the low-temperature regime. In contrast, it turns out 
that the photon subtractions before and after the noisy 
channel (thin solid and dashed curves) do not make any 
appreciable difference in the output entanglement. Actu- 
ally, two different strategies yield exactly the same states 
in the case of vacuum noise (nth = 0) and subtraction 
operation (i = 1, r = 0). In the case of low temperature. 




FIG. 3. Logarithmic negativity as a function of 77 = e"""', 
calculated in two different ways, estimation with covariance 
matrix (curves) and numerical calculation (symbols: circle, 
square, diamond), for three different states: no distillation 
(black thin curve and circle), coherent operation before noisy 
channel (blue thick curve and square), and coherent operation 
after noisy channel (red dashed curve and diamond). Squeez- 
ing parameter and thermal photon number are s = 0.029 and 
nth = 10~^, i.e., the cases of Fig. [2] (b). 



there is difference between two states, however negligible. 

On the other hand, when the thermal photon number 
of the reservoir increases (right panels), the operation 
before the noisy channel (solid curves) makes the output 
entanglement higher. This is closely related to the fact 
that a non-Gaussian entanglement can be more robust 
than Gaussian entanglement under a high-temperature 
reservoir, which was shown in Ref. [17| . In particular, 
it is crucial that entanglement can survive, even though 
very weak in its strength, through a noisy channel: Once 
entanglement dies out, there is no way of recovering it by 
LOCC. 

When we apply the distillation operation before a noisy 
channel, entanglement survives longer, e.g. with a small 
77 ^ 0.25 in Figs. [l](b) and (d). However, when we apply 
the coherent operation after the noisy channel, entangle- 
ment can survive shorter, e.g. 77 ^ 0.78 in (b) and rj ^ 0.5 
in (d). The separation time of TMSV, at which entan- 
glement completely disappears under a thermal reservoir, 
can be calculated by Simon's criterion [23|. Its analytic 
expression is given in |14l | by 



f-sep 



1 / 1 - e-2- 



(15) 



As nth increases, the separation time of TMSV becomes 
shorter, that is, Gaussian entanglement becomes more 
fragile under a noisy channel. If we perform distillation 
before the noisy channel, the resulting non-Gaussian en- 
tanglement can survive longer with a smaller rj. Non- 
Gaussian operations not only increase the degree of en- 
tanglement but also make entangled states more robust 
under noisy channels. We also note that the use of the 
coherent operation (thick curves) significantly enhances 
the performance in distillation over the photon subtrac- 
tions (thin curves) . As the squeezing of the initial state 
increases, however, we see that the difference in the dis- 
tillation performance between the two strategies becomes 
smaller [Fig. [5] (f )] . 



In order to better understand the advantage of non- 
Gaussian operation applied before the noisy channel, we 
analyze the logarithmic negativity in more detail for the 
states appearing in Fig. 2 (b) as examples. In particu- 
lar, we show in Fig. [3] both the degree of Gaussian-type 
entanglement (curves) that is obtained by the covariance 
matrix of a given state [Eq. (10)] and the degree of total 
negativity (symbols) for each state. 

Gaussian entanglement is very robust under a vacuum 
noise (or a very weak thermal noise) so that it can sur- 
vive the noise channel very long. Therefore it is deemed 
a good strategy for the purpose of entanglement distilla- 
tion that one applies a probabilistic operation after the 
noise channel. On the other hand, when the strength of 
thermal noise becomes rather significant, the Gaussian 
entanglement dies out at a certain value of interaction 
time [Eq. (15)]. For example, the initial Gaussian en- 
tanglement disappears at r/ ~ 0.8 in Fig. 2 (b). In this 
case, whatever operation is made on the evolved state 
after the noisy channel, the entanglement distillation be- 
comes impossible for rj < 0.8: Once entanglement dies 
out, distillation is impossible. Therefore one must per- 
form a certain operation before the state undergoes the 
noisy channel. 

As can be seen from Fig. 3, the non-Gaussian opera- 
tion (coherent operation) applied before the noise channel 
{rj = 1) increases both the Gaussian- type entanglement 
(arrow 2) and the total entanglement (arrow 3) from the 
initial Gaussian entanglement (arrow 1). With the in- 
creased entanglement due to the operation, the state 
becomes more robust than without operation. For ex- 
ample, the enhanced Gaussian-type entanglement (blue 
thick curve) alone evolves more robustly than the ini- 
tial two-mode squeezed state (black thin curve) because 
it has more entanglement (energy) that can resist the 
noise. Furthermore, due to the non-Gaussian operation 
performed on the initial state, there also exists a non- 
Gaussian type entanglement (roughly speaking, the dif- 
ference between squares and thick solid curve) which can 
also contribute to the surviving total entanglement. Note 
that even after the Gaussian-type entanglement disap- 
pears at 77 ~ 0.35, the non-Gaussian entanglement sur- 
vives until rj ^ 0.3. 

In a nutshell, at a significant level of thermal noise, the 
effectiveness of non-Gaussian operation before the noisy 
channel may be attributed to both the enhanced level of 
total entanglement (Gaussian and non-Gaussian type) by 
the operation and the robustness of non-Gaussian entan- 
glement in a certain parameter regime. However, such 
an effect becomes less remarkable when the initial entan- 
glement is rather high like the case of Fig. 2 (f). 



C. Success Probability 

From a practical perspective, it is important to con- 
sider the success rate for distillation operation as well as 
the degree of output entanglement. The success probabil- 
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FIG. 4. Distillation probability (in log scale) as a function of 
T] = e^""'. Legends and parameters are the same as those in 
Fig. [21 The curves are missing in the region where distillation 
is not possible at all, that is, negativity is before operations. 



ity of the coherent operation is equal to the probability of 
detecting a trigger photon, after erasing the which-path 
information on whether it is from the photon subtraction 
or addition, for each state. (See the proposed experimen- 
tal scheme in 20].) The success probability is shown for 
each case in Fig. S) We see that the distillation with a co- 
herent operation makes both of the ouput entanglement 
and the success probability higher than a mere photon 
subtraction. The success rate overall becomes an order of 
magnitude higher as the initial squeezing increases from 
s = 0.029, s = 0.114 to s = 0.403. 

In the low-temperature regime, as shown before, we 
can obtain the output state with high negativity by ap- 
plying a coherent superposition operation after sending 
the state through even a very lossy channel (dashed 
curves), e.g. ry ^ 0.1. However, the success probability 
is very low (< 10~^ ^ lO^'*) for 77 = 0.1 and it gener- 
ally increases with rj. On the other hand, if we apply 
the coherent operation before the noisy channel (solid 
curves), we have almost a flat success probability, ^ 10~^ 
in (a),~ lO^'^'* in (c), and ^ 10^^-^ in (e) regardless of 
V- 

In the high-temperature regime, the success probabil- 
ity for distillation after the noisy channel is higher than 
that for distillation before the channel. Especially, at 
77 ^ 0.7 in (d), both negativity and probability are higher 
when we apply operation after sending the state through 
the channel. However, for small rj where distillation after 
the channel is impossible due to the vanishing negativ- 
ity, we must apply the operations before sending the state 
through the channel. 



Note that the coherent operation is optimized for each 
case in order to maximize the degree of output entan- 
glement, not the success probabihty. The success proba- 
bihty is pfotted just according to the optimized coherent 
operation obtained that way, which can affect the be- 
havior of the success probabihty with respect to rj. For 
the case of operation apphed before the noise (blue solid 
curves) , the success probability by its definition is deter- 
mined before the state evolves under the noise. Thus, 
the case of photon subtraction (blue thin curves) must 
have a flat distribution with r]. In cases of Fig. 4 (e) 
and (f), where the degree of initial squeezing is rather 
high, it turns out that the optimized coherent operation 
is also very close to the photon subtraction regardless of 
T], thus it also shows a fiat distribution. On the other 
hand, for the case of operation applied after the noise 
(red dashed curves), the optimized operation varies with 
rj. For a rather high (small) 77, the optimized operation 
approaches the photon addition (subtraction) and the 
change of the success probability with 77 actually looks 
more prominent in linear scale than in log scale shown in 
the figures. 



IV. TELEPORTATION FIDELITY 
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FIG. 5. Teleportation fidelity as a function of transmitivity rj. 
Legends and parameters are the same as those in Fig. [2] The 
gray dashed line at F = 0.5 represents the maximal fidelity 
with separable resources. 



The logarithmic negativity studied in the previous sec- 
tion is a measure to quantitatively characterize the de- 
gree of entanglement. On the other hand, one may won- 
der how useful the distributed entangled state can be for 
informational tasks. We here investigate the usefulness 
of output entangled state by looking into the CV quan- 
tum teleportation [25|. Specifically, we investigate the 
teleportation fidelity, an operational measure of entangle- 
ment, to show how faithfully a given entangled resource 
can accomplish the quantum teleportation [25|. 

It is known that the fidelity between input and output 
states, averaged over all input coherent states, cannot 
exceed 1/2 without entangled resource [2a,[23|. With an 
entangled resource whose characteristic function is given 
by x(Cii^2), the input-output relation can be written as 
Xout(0 = Xin(e)x(r,0 m- The fidelity is then given 
by 



^=^/d'eXout(OXin(-0 



d'ex(r,Oxin(e)xin(-e) 



= ^/d'Cx(r,Oe- 



\i\' 



(16) 



-IClV2e««*- 



a*{ 



for a coherent-state input Xin(C) = e~ 

Using Eq. (15), we numerically calculate the telepor- 
tation fidelity as a function of 77 = e^'"* under two dif- 
ferent strategies, which is shown in Fig. [Sj Overall, the 
trend of the fidelity is very similar to that of the loga- 
rithmic negativity shown in Fig. [5] The teleportation 
fidelity can be high in the parameter region where the 



negativity is high, but this is not always the case. For 
example, in Fig. EJd), we find a crossover between two 
curves that represent the coherent operation before (solid 
curve) and after (dashed curve) the noisy channel, while 
there is no crossover in Fig. El^d). This implies that 
higher entanglement (negativity) does not always pro- 
vide a more faithful teleportation [2a]. From Fig. [51 we 
have a rather clear-cut conclusion about which strategy 
can yield a better performance for quantum teleporta- 
tion. In the low (high)-temperature regime, the coher- 
ent operation must be applied after (before) sending the 
TMSV through a thermal reservoir to obtain a higher 
teleportation fidelity. 



CONCLUSION AND REMARKS 



We have studied the entanglement distillation for CV 



states by a coherent superposition operation ta 



it 

when the TMSV decoheres under a thermal noisy chan- 
nel. We showed that a coherent operation still provides 
a greater advantage than a photon subtraction method 
in the presence of noise. In particular, we considered two 
different strategies for distillation: a coherent operation 
applied before and after the noisy channel. We investi- 
gated two entanglement characteristics, logarithmic neg- 
ativity and teleportation fidelity, to compare the perfor- 
mances of two strategies. In the low-temperature regime, 
where separation time of TMSV is rather long, the en- 
tanglement and the teleportation fidelity for the output 
state can be higher if we apply the operation after the 



noisy channel. The success probabihty for distiUation de- 
creases with the channel transmissivity 77 — e~^* . On the 
other hand, in the high-temperature regime, the entan- 
glement of TMSV decays rapidly and no distillable en- 
tanglement remains under a low transmissivity. It turns 
out to be better to apply operation before sending the 
TMSV through the channel so that the state becomes 
non-Gaussian with higher entanglement that can be more 
robust under a high-temperature noisy channel. In the 
limit of a large initial squeezing of the TMSV state, both 
strategies give almost the same degree of performance in 
distillation. 

In a realistic situation, the thermal photon number in 
environment can be very small. For example, the aver- 
age thermal photon for the wavelength A = 1064nm is 
TT-th — 2.61 X 10^^*^ at room temperature. In this case, 
the behavior of entanglement is similar to what we have 
shown with nth — 10^^. However, the average thermal 
photon number becomes as large as nth — 10^^ in the 
infrared regime (A = 20/xm) at room temperature and it 
becomes much larger in the microwave regime. We may 



also consider other kinds of noisy attenuator environment 
which can be modeled as a beamsplitter interaction with 
thermal photons. Therefore, a large photon number like 
nth — 10~^ may be considered as realistic. 

Finally, let us briefly address the applicability of dis- 
tillation with a coherent superposition operation for 
a long-distance quantum communication. In the low- 
temperature regime, we can obtain a highly entangled 
state when the coherent operation is successfully applied 
after the noisy channel. Thus the coherent operation 
alone can make a good distillation protocol to distribute 
entanglement between two distant nodes. However, in 
the high-temperature regime, transmission distance is 
limited and the degree of distributed entanglement is low. 
For a long-distance quantum communication, we may 
need to combine the coherent operation with other proto- 
cols such as quantum repeater and multiple-copy distilla- 
tion. Recently, CV distillation protocols have been stud- 
ied with single photon subtraction and two-copy mixing 
[2 3 . |30| . It may be worth further studying whether com- 
bining the coherent operation and two-copy mixing can 
be an efhcient distillation tool. 
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